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O ' 

O^l ' We analyze the behavior of unstable modes in the glasma produced in high energy heavy 

, ion collisions, using a simple model with effective homogeneous longitudinal color electric 

■ and magnetic fields. The unstable modes are approximately described as Nielsen-Olesen 

rj^ unstable modes under the homogeneous longitudinal gauge fields. We find that the Nielsen- 

Olesen unstable modes show properties very similar to those of the exponentially increasing 
unstable modes in the glasma recently shown by Romatschke and Venugopalan. Although 
initial gauge fields in the glasma are much stronger than the ones in our model, they decay 
with the production of the Nielsen-Olesen unstable modes. We discuss why we can reproduce 
the features of the glasma by using the effective homogeneous weak magnetic fields. Our 
analysis supports an idea that the decay of the gauge fields in the glasma is caused by 
Nielsen-Olesen instability. 
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§1. Introduction 

In high energy collisions of heavy ions in RHIC or LHC, the most important 
ingredients for producing quark gluon plasma ( QGP ) are small x gluons in the 
nuclei. x )~ 5 ) The small x gluons with also small transverse momenta are sufficiently 
dense in the nuclei and so they may be treated as classical fields produced by glassy 
large x gluons even after the collisions of the nuclei. It has been shown that longitu- 
dinal color magnetic and electric fields of the small x gluons are generated initially at 
the collisions. They are classical fields and evolve classically according to a color glass 
condensate (CGC) model. 6 ) These fields, which are called as glasma, are expected 
to give rise to QGP by their rapid decay. 7 )' 8 ) 

In an extremely high energy collision, the thickness of nuclei is nearly zero due 
to the Lorentz contraction. Thus, the initial gauge fields have only transverse mo- 
mentum perpendicular to the collision axis, but have no longitudinal momentum ( 
rapidity dependence ). Such classical gauge fields can not possess any longitudinal 
momentum in their classical evolution, since equations of motion of the gauge fields 
are invariant under the Lorentz boost along the collision axis. 2 ) -4 ) Thus, the naive 
application of the CGC model, e.g. McLerran- Venugopalan model ( MV model ), 6 ) 
does not give rise to thermalized QGP. It has recently been shown, 9 ) however, that 
the addition of small fluctuations with rapidity, e.g. quantum fluctuations 10 )' 11 ) 
to the initial gauge field induces exponentially increasing modes with longitudinal 
momentum. The production of the exponentially increasing modes implies that a 
process toward thermalization has started; the decay of the gauge field and the 
isotropization of momenta. Although the decay of the classical fields have been 
clarified in the numerical calculations, the physical mechanism of the decay is still 
unclear. 
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In this paper we show using a simple model of an initial gauge field that the 
decay of the glasma is caused by Nielsen-Olesen unstable modes. 8 )' 12 )~ 14 ) The modes 
arise under the initial gauge field when small fluctuations around the gauge fields are 
taken into account. Especially we compare the time evolutions of the unstable modes 
with those of the exponentially increasing modes of longitudinal pressure shown in 
the previous calculation. 9 ) We find that the previous main results can be reproduced. 

The initial gauge field in our model is spatially homogeneous and much weaker 
than the fields of the glasma. 9 ) The glasma are strong ( e.g. color magnetic field 
B ~ 0(Q 2 S / g) ) and inhomogeneous ( their coherent length being the order of ~ Qj 1 ; 
Q s is a saturation momentum written as Q s = g 2 fj, in the ref. 9 ) ) The initial gauge 
field is chosen to reproduce the time evolutions of the instabilities observed in the 
reference. 9 ^ It can not reproduce local behaviors of the instabilities in the transverse 
space since it is homogeneous in the space. 

Nielsen-Olesen instability is associated with color magnetic fields, not color elec- 
tric fields. In the paper we are only concerned with the instability of the color 
magnetic fields. The color electric fields produced in heavy ion collisions are unsta- 
ble against quark pair creations 15 ) so that they may decay sufficiently fast. We do 
not discuss the instability of the color electric fields. 

In the next section, we explain the instabilities observed in the previous sim- 
ulation. 9 ) In the section(3), we review briefly Nielsen-Olesen instability. In the 
section(4), we explain the assumption used in our model, expecially the relevance of 
the use of homogeneous color magnetic field instead of inhomogeneous glasma. In 
the section(5), we introduce our simple model for analyzing the decay of the glasma 
in r and rj coordinates. In the section(6), we compare our results with those obtained 
in the previous simulation. 9 ) In the final section, we conclude our results. 

§2. Instabilities in initial gauge fields 

In high energy heavy ion collisions, classical color gauge fields are generated and 
are pointed to the collision axis ( longitudinal direction ). They are color electric 
and magnetic fields. They arise at r = just when the collisions occur. ( Here we 
assume collisions in high energy limit where heavy ions are Lorentz contracted to 
have zero width. Thus, the collisions occur at the instance of r = 0; r = \J x\ — x\. 
) They have sufficiently large energy densities to produce thermalized quark gluon 
plasma in the subsequent their decay. The color gauge fields are coherent states 
of small x gluons with transverse momenta less than a saturation momentum, Q s . 
These gluons are described by a model of color glass condensate, e.g. MV model. 
More explicitly, the color gauge fields are given initially at r = as functions of 
color charge density of large x gluons inside of heavy ions; the color charge density 
is determined with a Gaussian distribution in the MV model. 

Such classical gauge fields are uniform in the longitudinal direction. On the 
other hand, they are not uniform in the transverse directions; the scales of their 
variations in the directions are typically determined by the saturation momentum 
Q s . This is because they are made of the gluons possessing transverse momenta 
typically given by Q s . This implies that the fields keep their directions ( parallel 
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or anti-parallel to the collision axis ) inside of transverse regions whose widths are 
given typically by Qj 1 . But, they change their sign outside of the region. In this 
way their directions are never uniform in the transverse directions, although they 
are uniform in the longitudinal direction. 

Anyway, such color gauge fields are given initially in the collisions. After their 
production they evolve according to the gauge field equations. Since the gauge fields 
equations are invariant under Lorentzs boost along the collision axis denoted by xz 
( rj = log ( j 1j| ) —> V + cont. ), the color gauge fields have no dependence of rj in 
their development. That is, the gauge fields are still uniform in the longitudinal 
direction after the collision. In such a circumstance, a numerical calculation 5 ) has 
been performed to show that they evolve smoothly in time r and become weak 
owing to their expansion. Any unstable behaviors of the fields were not found in the 
calculation. Non-linearity in the gauge field equations seems not to play dramatici 
roles in the evolution of the gauge fields. Actually, it has been shown 14 ) that their 
time developments can be reproduced qualitatively in the linear analysis of the field 
equations; especially the fields become weak with time as 1/r. (A kind of "Abelian 
dominance" 16 ) seems to hold for large r since self-interactions of gauge fields become 
ineffective because of smallness of the fields in the large r. ) We use the feature in 
order to make our simple analytical model of the glasma decay. 

Subsequently, a numerical simulation has been performed 9 -* with a slightly mod- 
ified initial condition. That is, much small gauge fields depending on r\ are added by 
hand to the original initial gauge fields. Such gauge fields may arise because the uni- 
formness in r\ of the initial condition is broken in real situations: Heavy ion collisions 
occur with finite energies or the initial conditions derived in the MV model receive 
high order quantum corrections. Both of them give rise to much small corrections 
depending on rj to the initial gauge fields without rj dependence. 

The simulation has clarified the existence of unstable modes in the evolution of 
the gauge fields; Fourier components in 77 of longitudinal pressure increase exponen- 
tially in r. This implies that a component of gauge fields added in the simulation 
increases exponentially. The existence of the modes implies that the initial color 
electric and magnetic fields uniform in the longitudinal direction are unstable under 
the small fluctuations depending on rj. 

We should note that since gauge field fluctuations depending on rj are suffi- 
ciently smaller than the initial gauge fields, they can be treated perturbatively. 
That is, the fluctuations evolve under the background initial gauge fields without 
self- interactions. Hence, the analysis of their evolution can be performed in the 
linear approximation. 

Here, we review characteristic properties of the unstable modes found in the 
numerical simulation. We denote the longitudinal pressure as P^. When we denote 
Fourier components as P v (k v ,T) = J drj P v (rj,r) exp(ik v rj) ( k v denotes longitudinal 
momentum ), P v (k v ,T) evolves smoothly in r when initial gauge fields have no de- 
pendence of the rapidity rj. Obviously, P^(fc^,r) oc 5(k v ) in the case. Once gauge 
fields depending on rj are added to the initial gauge fields, P r) (k ri ,T) shows exponen- 
tial increase. Namely some unstable modes are excited when gauge fields fluctuations 
depending on rj are added. P v (k v , r) shows the exponential increase just after a time 
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r{k rj ) > passes. In other words, P v (k v ,T) does not increase exponentially until 
the time r = T(k v ). The simulation shows that r(k v ) oc k v . That is, the compo- 
nent Pti(t, k„) begins to increase exponentially later than the time when components 
Pt](t, k' v ) ( k' < krj ) increase. Hence, we can define a maximum momentum ^(max) 
at the instance r such that the component P v (t, fc^(max)) begins to increase expo- 
nentially at the instance; components P v (t, k v < fc^(max)) have already increased 
exponentially. The simulation shows that £1,7 (max) increases linearly with time r. 

Another interesting property found in the simulation is concerned with longi- 
tudinal momentum distribution at a time r, namely P r) (k r) ,r). The distribution 
has a peak at the value of k v (p), which is almost independent of r. This implies 
that the unstable modes with the characteristic longitudinal momentum k v (p) are 
induced dominantly. The peak in the distribution, P v (k v (p),T), increases such as 
P v (k v (p) , t) oc exp(const.r 1 / 2 ). 

We should mention that the growth rate of the peak is much smaller than the 
typical time scale Qj 1 in the system. Thus, it takes a too long time for the small 
gauge field fluctuations to become comparable order of magnitude of the initial gauge 
fields. This is a bad new for the realization of thermalized QGP through the decay 
of the gauge fields. 

We will show that these characteristic properties can be understood using Nielsen- 
Olesen instability. The point in our analysis is that the background initial gauge fields 
are approximately described by Abelian gauge fields and the gauge field fluctuations 
added can be treated perturbatively. Based on the approximation, we show that the 
unstable modes found in the numerical simulation are just Nielsen-Olesen unstable 
modes. 

§3. Nielsen-Olesen Instability 

We review briefly Nielsen-Olesen instability by using SU(2) gauge theory. The 
instability means that homogeneous color magnetic field is unstable in the gauge 
theory. In order to explain it, we decompose the gluon's Lagrangian with the use 
of the variables, "electromagnetic field" = A^, and "charged vector field" = 
(A^ + where indices 1 ~ 3 denote color components, 

L = -\$%, = -\{d»A v - d v A»f - \\D^ V - D v ^\ 2 - 

+ ieid^ - d v A^\$ v + ^(<?V<4 - <Pu^l? (1) 

with = <9 M + igA^, where we have omitted a gauge term = 0. 

We can see that the charged vector fields ^ couple with electromagnetic field 
A^ minimally through the covariant derivative and non-minimally through the 
interaction term, ig{d^A v — d u A fl )<pl l <P u . When a homogeneous color magnetic field 
B > described by = A^ is present, we analyze the fluctuations 0^ under 
the color magnetic field. To do so, we solve the following equation of under the 
background field B, 
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(d 2 -D 2 T 2gB)<l>± = (2) 

with Af = (-Bx 2 ,Bx 1 ,0)/2, where c/>± = (<Pi ±i<P 2 )/V2 and D = d + igA B . Index 

± denotes a spin component parallel (+) and anti-parallel (— ) to B = (0,0, B). 
We have assumed the magnetic field pointed into 23 direction. In the equation, 
higher order interactions have been neglected. The energy spectra of the fields cp± 
are easily obtained. The energy u of the charged vector field <p± oc e lu)t is given by 
J 2 = k 2 + 2gB(n + 1/2) ± 2g.B. Integer n > denote Landau levels and does 
momentum parallel to the magnetic field. 

The term ±2gB in lo 2 comes from the non-minimal interaction which represents 
anomalous magnetic moment of the charged vector fields. It is obvious that the 
modes of 0+ have imaginary frequencies oj 2 = k 2 — gB < when n = and k 2 < gB. 
The modes occupy the lowest Landau level ( n = ). This implies that the field 0+ 
increase exponentially in time. The modes with ui 2 < are called as Nielsen-Olesen 
unstable modes. Therefore, when homogeneous color magnetic fields are present, the 
states are unstable; the Nielsen-Olesen unstable modes are generated spontaneously 
and the states decay into more stable states. This is the Nielsen-Olesen instability 
in the gauge theory. 

For the convenience of later discussions, we write down the Hamiltonian 17 ) of 
the field (/)+ neglecting higher order interactions, 

H = \d t cf> + \ 2 + \Dcf )+ \ 2 -2gB\ ( j) + \ 2 , (3) 

where the term |_D(^>_|_| 2 becomes (k 2 + gB)\cj) + \ 2 when cf) + occupies the lowest Landau 
level under the homogeneous color magnetic field B. The Hamiltonian holds even 
for inhomogeneous magnetic field B = eijdiAj. We can see that even when inhomo- 
geneous magnetic field B is present, the presence of the field 0+ can make lower the 
energy than the energy ( = )of the state with 4> + = 0. In the sence, inhomogeneous 
magnetic fields are also unstable in general. ( In order to demonstrate the instability 
we need to show the presence of bound state solutions by solving eq([l]) in the next 
section. ) We call the field (/>_)_ as Nielsen-Olesen field. Similarly, the field 4>- is a 
Nielsen-Olesen field for the inhomogeneous B since it also describes unstable modes; 
the Hamiltonian of </>_ is given by H{<f)_) = \dt<f)+\ 2 + |i5(/>+| 2 + 2gB\4> + \ 2 and the 
sign of gB can be negative in the transverse space. 

§4. Our assumption for analysis of instability in gauge fields 

First of all, we explain assumptions used in our simple model for gauge field 
evolution in heavy ion collisions. Here we use the Cartesian coordinate for the 
explanation. Initial gauge fields independent of the rapidity used in our model are 
only longitudinal color electric and magnetic fields. They are maximal Abelian 
components of gauge fields. For example, Af = (0,0,^4?) in the SU(2) gauge theory. 
Then, non- Abelian interactions vanish and only linear equations remain. It has been 
discussed 14 ^ that even if we make such a simplification of initial gauge fields in the 
glasma, we can reproduce quite well the numerical results 5 ^ in their evolution at least 
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for r > Qj . Thus, the approximate use of such Abelian initial gauge fields instead 
of the non-Abelian glasma is appropriate. This is the first our assumption. 

Such Abelian gauge fields are inhomogeneous in transverse space similarly as 
the glasma. Then, as we have discussed in eq([3]), the color magnetic fields B are 
unstable energetically owing to the presence of Nielsen-Olesen fields. This instability 
represents instability of the glasma in our model. It is described by the equation, 

J 1 ^ = (_L>2 + jfejj - 2gB)(j) (4) 

with Dj, = (dT+igAx) 2 and B = rot At, where we assume that <f> oc exp(iu;t — 2/03X3). 
The equation can be derived from the Hamiltonian in eq©. It is a "Schrodinger 
equation" of charged particles under the magnetic field B with an additional potential 
term —2gB. There are solutions of "bound states". The binding energy is given by 
—uj 2 > 0. Thus, the growth rate of the unstable modes is given by the imaginary 
part of the frequency u. 

The second our approximation is to replace inhomogeneous B in the operator, 
— D T — 2gB in eq(j3]) with effective homogeneous one B, keeping the eigenvalue of 
the operator —D T — 2gB unchanged. We should note that the dependence of the 
longitudinal momentum in the growth rate and the momentum distribution ^(k^) 
do not change even if we make the replacement as far as eigenvalues of the operator, 
— D T — 2gB, are the same. 

Thus, by the replacement we can analyze what we are concerned with, that 
is, the time evolutions of the unstable modes and the characteristic features, e.g. 
growth rates, associated with the evolutions. But local behaviors of the fields in the 
transverse space are lost by the replacement. The replacement is possible in principle, 
but is difficult in practice. In our simple model, the value of B is determined by 
comparing our results with the previous simulation. 9 ) 

§5. Our simple model of Nielsen-Olesen instability in glasma 

Now we explain the detail of our simple mode for the glasma decay. We use the 
proper time coordinate, r = \Jx\ — x\ and rapidity, 77 = log( ^°^| ), as a longitudinal 
coordinate. The coordinate is convenient for the description of expanding "glasma" 
generated in the early stage of heavy ion collisions. The collisions occur at X3 = ( 
rj = ) and xq = ( r = ) in extremely high energies. Thus, the heavy ions are 
Lorentz contracted to have vanishing width in the longitudinal direction and only 
extended in the transverse directions with coordinates Xi = (x\,X2)- 

We analyze SU(2) gauge fields, A^. Corresponding gauge fields in the coordi- 
nates are given such as A T , A v and Ai = (Ai, A2). Taking a gauge condition, A T = 0, 
the Lagrangian of the fields is given by 

tL = r (^{d T A v f + ^(drAf - - Ip^j , (5) 

with F 2 - = (drjli - diA v + gA v x Ai) 2 and Ffa = (diAj - djAi + gli x Aj) 2 . 
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We define the complex fields, (ft v , <pi, and real ones, A v , A{ by rearranging color 
components of the gauge fields, 

Ar t = {A^, A 2 , A^) = ( v v , v ^^ / - V ,A ri ) and A4 = ( — * , i , Ai), 

(6) 

and define the following "charged field" with spin parallel, (ft + ( anti-parallel, ) to 
the longitudinal direction, </>± = ^ 1 ^' 2 ■ Initial gauge fields are introduced as maxi- 
mal Abelian component A v , Ai of the gauge fields as we mentioned above. It is easy 
to see that the fields, 4>r) and <p±, are transformed under the gauge transformation, 
A^ — > WApU + g~ 1 Wd fl U with £7 = exp(i#(73), such that (f) Vj ± — > exp(— iO)(f> Vj ±, 
where <7j are Pauli metrices and 6* is a constant independent of r. Therefore, we 
may think that the complex fields, 4> v ^± are U(l) charged fields corresponding to the 
symmetry. 

We introduce a longitudinal color magnetic field B and an color electric field E 
as initial background gauge fields, both of which are assumed to point to the third 
direction in SU(2) gauge group; B = eijdiAj and E = ^drA^. The fields point 
to the direction parallel to the collision axis in the real space; B = (0, 0, B) and 
E = (0,0, E). We assume that the background fields are generated at r = in 
heavy ion collisions. 

With the use of the fields, Ai ;q and (f)± :V , the Lagangian leads to 

rL = T -(d T Mf + ±(d T A v ) 2 + r(|^ + | 2 + \dA-?) + l -\dr^\ 2 ~ \fh ~ ±fli 

- r(|A0+| 2 + IA0-I 2 ) - ^(lA^+l 2 + |A^-| 2 + |A^| 2 ) + 2t 9 B(\^\ 2 - |0_| 2 ) 
+ T -\D_<j) + + D + ^\ 2 + -=-((£>_0+ + D+^D^r, + c.c.) + 

with /jj = djA,- - <9jA, /^j = - SjA/, = fr h i - if v ,2, A = + 

D v = d v + igA v , and At = A ± *A, where we have neglected surface terms just 

like djJj. 

Obviously, the Lagrangian is invariant under the U(l) gauge transformation, 
0±,r? —> 0±,r? exp(— z#) along with — > + g~ 1 di tV 6. The kinetic energies of the 
fields Ai, A v , <p± and are presented in the first line of the Lagrangian. In the 
second line, minimal interactions between <f>±^ and the abelian gauge fields, Ai and 
A v are presented. We can see in the second line that the charged fields <f>± receive 
anomalous magnetic moments, 2rgB(\(p + \ 2 — \4>-\ 2 ). This term plays an important 
role for making the field 0+ unstable, that is, Nielsen-Olesen unstable mode. There 
are the quartic interactions of the fields in the fourth line, which describe repulsive 
forces among the fields <fi v ,±- The repulsive force leads to the saturation of the 
exponential increase observed in the simulation. 9 ) The terms in the third line are 
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irrelevant to our discussion below ( the terms can be gauged away ). 

When an initial gauge field configuration of B = eijdiAj and E = ^d T A v is 
given at r = 0, the subsequent evolution of the fields A v and Ai is governed by the 
equations, 

d T {^d T A v ) - ^(dfA v - dr^Ai) = o, 

and d T (Td T Ai) - r(d]Ai - OidjAj) + (<9 2 ^ - did v A v ) = 0. (8) 

The equations are obtained from the Lagrangian by neglecting the other charged 
fields. The approximation is valid when the other fields are sufficiently small so 
that the interactions between E, B and the others can be neglected. A solution 
of B and A v independent of the rapidity rj is given by B = BqJq{Qqt) cos(Qqx) 
and A rj = ct,J\{Qqt) cos(Qox), where we suppose that the fields carry a transverse 
momentum, Qq ( Qq = \Qq\ ) as a typical transverse momentum of backgound gauge 
fields. Bq is a constant and Jq,\{Qqt) are Bessel functions. ( General solutions are 
given by the average in Qq over momentum distritutions. ) A constant of c may 
be determined by the requirement that B = E as r — > 0. It leads to c = Bq/Qq. 
The requirement arises from the initial condition of (Tr(£? 2 )) = (Tr(E 2 )) for r — > 
given in the MV model of CGC. 18 ^ Here, the expectation of (~) is taken over the 
distribution of large x gluons according to the MV model. 

Here we neglect spatial dependence of the gauge fields according to the second 
our assumption. Additionally we simplify the factor Jq(Qqt) such that Jq{Qqt) — ► 
sin(Qo' r ) \/2/kQqt ~ \/^/^Qqt by neglecting the oscillating factor sin(Qo7"). This 
smooth decay of the fields roughly coincides with the numerical evolutions of the 
glasma. 5 ) Hence, the simplification is appropriate for discussing small fluctuations 
around the slowly decaying background gauge fields. Therefore, we assume the 
following background initial gauge fields, 

S = SoV2/vtQ t and A v = B (t/Q ) ^2/ttQ t, (9) 
which reproduce the smooth decay of (Ti(B 2 )) and (Tr(.E 2 )) for large r. 

Under the background gauge fields, we analyze the development of the small 
fluctuations 4> v ,±- The fluctuations correspond to the small fluctuations added to 
initial background gauge fields in the previous simulation. 9 ) Since they are supposed 
to be much small, we take into account only quadratic terms of the fields in the 
Lagrangian. These fluctuations in general oscillate with small amplitudes. But one 
of these fluctuations increases exponentially. In our model the field (f) + ( or ^_ when 
gB < ) is the one increasing exponentially with r. The other fields simply oscillate 
and their amplitudes remain small. 

We write down the equation of motion of the field 4> + , 

d 2 T <t> + + \d T ct> + + ^~9Mr)f _ gB{r) ^ ^ = Qj (1Q) 

where k v denotes a longitudinal momentum; + oc ex.p(—ik v ri). We have taken only 
a component in the lowest Landau level. It is easy to see that due to the last term 
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—gB(r), 0+ increases exponentially just as </>+ oc exp(\/gBr) for r — > oo when gi? 
is independent of r, i.e. a solution of the equation, <9 2 </>+ — gB(f>+ ~ 0. Their wave 
functions are given such that (fi + = g m (T)z m exp(— \x\ 2 /4:l B ) with z = x% + ix2 and 
integers, m > 0, where Ib = 1/y/gB denotes cyclotron radius and flVnfr) is governed 
by the eq (|10p ; we have neglected the smooth expansion of the cyclotron radius. 

§6. Our results 

In order to solve the eq (|10p with A v (t) and B{t) given in eq©, we rewrite the 
equation in the following, 

dU + + + ^~*f ? - 7=7)*+ = °> C 11 ) 

where dimensionless parameters are defined as r' = Q s t, a = \j2j^i y gB^IQ\) x 
(Qo/Qs) 3 ^ 2 and b = y / 2/-K(gB /Q' 2 ) ) x (Qo/Qs) 1 ^ 2 ■ In the subsequent calculations 
we treat the scale of the field (/>+ arbitrary although it is much smaller than the 
background field. This is allowed in the approximation of taking only quadratic 
terms of the field 0+ in the Lagrangian. The coefficients a and b are determined 
for reproducing the results in the previous simulation; 9 ) actually we have used a = 
(0.05) 2 and b = 0.38 in order to obtain our curve in Fig. 5. 

Before solving the equation numerically, we briefly explain how the solutions 

behave with r'. The term of ^ ( fc '> i^lJ _^L, J ( = u 2 ) is just a spring constant. 

The term of \d T i(f>+ represents a friction, which becomes weaker as r' becomes 
larger. Thus, we understand that the field oscillates as far as w 2 > 0, that is, in the 
early stage ( r' ~ O(l) ) after the production of the background fields. The spring 
constant a; 2 becomes small with r'. Once w 2 becomes negative, the field stops the 
oscillation and begins to increase exponentially. 

In Fig.l we show the typical behavior of the field, <p+{k ri = 16.5, r'), with the 
initial conditions of </>+(t' = 0.01) = 1 and 9 T + (r' = 0.01) = 0. ( Obviously, 
taking different initial conditions do not change the global behavior of <j) + for large 
t', since it simply oscillates in small r'. ) The field increases exponentially after 
the oscillation in the early stage. When the longitudinal momentum k v becomes 
large, the time r when the field 0.j-(& 7 j,t / ) begins to increase exponentially, becomes 
large. This implies that as r' becomes larger, the modes with larger longitudinal 
momentums are excited. 

In Fig. 2 we show the time dependence of the maximal momentum L(max). 
The maximal momentum ^(max) at the time r' is defined as the momentum with 
which the mode <f>+ (k v (max), r') starts to increase exponentially at the time r'. The 
modes with fc„ < fc„(max) have already been increasing exponentially at the time 
t' . The maximal momentum ^(max) may be obtained by solving the condition 
of the spring constant u; 2 = 0, but /^(max) is defined as |0+(fc^(max), t')\ = 2 in 
our calculation. We have shown both results in Fig. 2, which almost coinside with 
each other, ^(max) increases almost linearly in r', but the solution of u> 2 = 



10 



Aiichi Iwazaki 



shows that ^(max) oc r' 3 / 4 . The result agrees with the previous one, 9 ) although 
the rate of the increase is approximately four times smaller than the previous one; 
fe^(max) ~ 0.015 Q s t + 5. 

It has been shown 9 ) that ^(max) deviates from the linear dependence around 
the time when the exponential increase is saturated. After the deviation, k v (m&x) 
increases very rapidly. These phenomena could be understood in our model as the 
result due to the onset of quartic interactions among the field <fr+. The quartic 
interactions make the energy of the mode with largest amplitude be transmitted to 
the other modes with higher longitudinal momenta. 




Before proceed to make the comparison, we should mention why we compare 
the evolution of the field ^+ with the evolution of longitudinal pressure discussed in 
the previous simulation. 9 ) Their behaviors are quite similar to each other. Roughly 
speaking, there are two components of fields involved in the ref., 9 ) large ones and 
much small ones. The large ones ^4(large) are just the background fields independent 
of rapidity, which are produced according to the MV model. On the other hand, the 
small ones a(small) are fluctuations depending on the rapidity which are added by 
hand to the background fields. In the circumstance, the longitudinal pressure is 
composed of two parts, P Vj o and 5P V ; P r] = P rh § + 5P n . 

P v = t~ 2 ($% ti + (rdrAif^j - F* 2 - (^d T A v ) 2 ~ P vfi + SP V , (12) 

where P Vt o is rapidity independent and is formed only of the large components, 
while SPrj is formed of the small components as well as the large ones; 5P V = 
P(A(large)) x a(small). The linear dependence on a(small) comes from the approxi- 
mation of a(small) <C A(large). P V fi decreases smoothly with r. On the other hand, 
5P V increases exponentially although it is still much smaller than -P^o- Therefore, 
some of small components a(small) depending on the rapidity, increase exponentially 
as has been shown. 9 ) Our simple model of gauge field evolution indicates that such 
a small component is just the Nielsen-Olesen unstable mode (/>+. That is the reason 
why we compare the evolution of the field (j) + with the evolution of the longitudinal 
pressure. It should be noted that the Fourier component of SP^ in the rapidity is 
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Fig. 3. distribution of longitudinal momen- 
tum k v at r' = Q s t = 1500. 
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Fig. 4. t' dependence of k v (p) ( dot ) and 
k n = 0.18^ ( solid ). 



determined only by the factor of a(small), which corresponds to small fluctuations 
<t>r),± in our model. 

We now proceed to make the further comparison. In Fig. 3, we show the typical 
momentum distribution ^-^{k^, r')\ at r' = 1500. The distribution in k^ has been 
obtained by solving eq (|lip for 4> + (k, 1 , r') with each k v chosen within the range 0.24 < 
k v < 18 by 5k v = 0.03. The distribution is not smooth but oscillating rapidly in k v . 
( When we magnify a small region, e.g. Sk^ ~ 0.2 of the Fig. 3, we can see explicitly 
the oscillation in k^. Roughly speaking, this oscillation comes from the oscillation 
of a modified Bessel function Ik^iQsT). ) But we can find out a smooth distribution 
obtained by averaging the original one over an appropriately small but sufficiently 
large Sk^ to be compared with wave length in the oscillation. Then, the smooth 
one coincides almost with the mountain-like form of the distribution in Fig. 3. ( The 
average corresponds to an average over initial conditions for solving eq (|ll|) . This is 
because slightly different initial conditions lead to slightly different curves. ) We 
find that the distribution has a peak at a momentum fc„(p). The mode with the 
momentum is generated most efficiently. 

In Fig. 4, we show how k v (p) depends on the time t ; k„(p) increases very slowly 
with t' just as O-ISV^- Furthermore, we can see that the momenta are smaller 
than /c,j(max). These results agree with the previous ones. 9 ) The presence of the 
longitudinal momentum ^(p) — 6 ~ 8 almost independent of time implies that 
in the decay of the gauge fields uniform in rj, specific modes with the momentum 
ferj(p) is generated most efficiently, which breaks the homogeneity in rj. We do not 
understand why such specific momenta are present. 

Finally, we show in Fig. 5 how \<j) + (k v (p), r)| increases with r', that is, the time 
dependence of the peak |</> + (fc^(p), r)|. It increases as exp(r' 3 / 4 ) for r' — > oo, which 
can be read from the equation dTTJ. The growth rate is defined as (log \ 4>\)/t. Thus, 
the growth rate of the field \4> + (k v (p), r)| decreases slowly with r. Since k v (p) is 
almost constant in time, \4>+(k v , r)| also increases in the similar way to 1 7r_|_ (A;^ (p) , -r) | 
for any k v . 

For a comparison, we have depicted the Fourier component of the longitudinal 
pressure, P v (k v (p)) = (do + d\ exp(0.427\/r')) in the ref. 9 ) as well as the function 
of exp(0.00544r') also used in the ref. 9 ) Obviously, the field, \4> + (k v (p), t')\, in 
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Fig. 5. longitudinal pressure P v (k v (p)) (solid 
) in ref. 9 ' , \4>+\ ( short dashing ) and 
exp(0.00544r') ( dashing ). 

our calculation agrees with the longitudinal pressure better than the function of 
exp(0.00544"r'). ( In Fig. 5, we take the scale arbitrary in the vertical coordinate. In 
order for \(f> + (k ri (p),T')\ to coincide precisely with the pressure shown in the ref., 9 ) we 
have only to take the scale of the field \(j) + (k v (p), t')\ appropriately. ) Our simulation 
does not reproduce strictly the behavior of the pressure such as exp(v / r')- We expect 
that the more elaborate treatment of the background magnetic field may give rise 
to the behavior exp(V / r')- ( When B decreases as ~ 1/r' instead of 1/\/t', this 
behavior can be obtained. We will discuss the validity of the behavior B ~ 1/r' in 
near future. ) 

As it has been shown, the behaviors of the longitudinal pressure calculated in the 
MV model with the small fluctuations added can be roughly reproduced in our simple 
model: 1) ^(max) increases linearly with r, 2) fc^p) depends on r very weakly and 
k v (p) is much smaller than k v (m&x), and 3) the pressure P v (k(p)) increases with r 
as exp(r 3 / 4 ), although the pressure obtained in the simulation increase as exp(y / r). 
Furthermore, we may argue that the saturation of the exponential increase arises due 
to the repulsive self-interaction of (fi+ in our model. The repulsive interaction need 
more energies for the field to increase more. Thus, it would stop the field increasing. 
Although our results are different in detail with those in the simulation, the rough 
agreement shows that our simple model of instabilities is valid approximation for the 
instabilities observed in the glasma. Thus, the instabilities of the glasma observed 
in the previous simulation is caused by the Nielsen-Olesen unstable mode. 

In order to obtain these results, we have used the parameters such as a = (0.05) 2 
and b = 0.38. Roughly speaking, the parameter o gives the growth rate of the field 
cf) + . Thus, the parameter can be determined by making it fit the growth rate obtained 
in the simulation. 9 - 1 But the growth rate shown in the simulation is the one of the 
field (p + (k ri (p)), not cp + (k ri ) itself. Here k v (p) depends on r although the dependence 
is very weak. This requires a careful adjustment of the other parameter b. 

The determination of the parameters in detail is not important. An important 
thing is that these parameters lead to a weak homogeneous color magnetic field B$. 
Actually, the above values of a and b corresponds to the physical parameters of 
Qo — Qs/152 and gBo ~ 4.7Qq. Thus, the growth rate ~ \fgB~Q discussed in the 
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section (4) is much smaller than y^glB] ~ Q s . Indeed, the growth rate is about 
0.00544(3 <j derived from the reference function exp(0.00544r') depicted in Fig. 5. 

We should mention that the smallness of the growth rate comes from the inho- 
mogeneity of the glasma. The potential —2gB ~ Q 2 for the field </>+ in eq§%$) varies 
so rapidly in transverse space; it has many attractive regions ( —2gB < ) and 
repulsive regions ( —2gB > ) whose widths are of the order Qj . Wave functions 
of the bound states are extended over these regions involving both attractive and 
repulsive potentials; they can never be trapped within a region with attractive po- 
tential. Therefore, the binding energies —oj 2 become much smaller than Q 2 . Thus, 
the growth rate becomes much smaller than Q s . 

§7. Conclusion 

To summarize, we have discussed that inhomogeneous color magnetic fields 
gB ~ Q 2 produced in the high energy heavy ion collisions decay with the pro- 
duction of the Nielsen-Olesen fields Instead of analyzing the evolution of the 
field 0+ under the inhomogeneous color magnetic fields, we analyze it using an effec- 
tive homogeneous color magnetic field. Then, we have compared the time evolution 
of the Nielsen-Olesen field (/>+ with the evolution of the longitudinal pressure shown 
in the simulation. 9 ) We have found that our simple model with the effective weak 
homogeneous color magnetic field gB$ reproduces the important features clarified 
in the simulation; growth rates, longitudinal momentum distributions, etc. of the 
unstable modes. The coincidence is not accidental. These are properties associated 
with time and longitudinal directions. Even if we use the homogeneous magnetic 
field, the properties of the unstable modes can be reproduced in principle. On the 
other hand, properties associated with the transverse directions can not be repro- 
duced with the use of such fields. Therefore, our analysis shows that the decay of the 
glasma generated initially at r = in heavy ion collisions is caused by Nielsen-Olesen 
instability. 

On the other hand, there is an insistence 9 )' 19 ) that Weibel instability known in 
plasma physics is the cause of the glasma instability: Inhomogeneous electromagnetic 
plasma whose momentum distribution depends only on transverse momentum, shows 
Weibel instability under small magnetic field applied. The instability is discussed 
by using the Boltzmann equation of charged ( color charged ) particles coupled with 
electromagnetic ( color gauge ) fields, while the glasma instability in the simulation 
has been shown in pure gauge theory. Nielsen-Olesen instability is the instability 
arising in pure gauge theory. In this sense, the relevance of the Weibel instability to 
the glamsa instability is not obvious. Thus, it is reasonable to think that the glasma 
instability shown in the simulation is just Nielsen-Olesen instability. We will discuss 
in future publications why the glamsa instability is just Nielsen-Olesen instability, 
not Weibel one. 

We would like to express thanks to Dr.K. Itakura in KEK and Dr.H. Fujii in 
University of Tokyo for useful comments. 
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